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ABSTRACT 

Context. Ultraluminous infrared galaxies are merging systems characterized for containing large amounts of dust and emitting the bulk 
of their energy in the infrared. Dust affects several of the observed properties derived from optical and near-IR rest-frame data such 
as the stellar morphology, and the ionized gas distribution. Systematic analysis of the dust distribution in representative samples of 
ULIRGs are needed to investigate its two-dimensional structure, and to establish its impact in the derivation of fundamental properties 
such as star formation rates, effective radii, and dynamical masses. 

Aims. We investigate the two-dimensional kpc-scale structure of the extinction in a representative sample of local ULIRGs using the 
Hof/H/? line ratio. 

Methods. We use optical integral field spectroscopy obtained with the INTEGRAL instrument at the William Herschel Telescope. 
Complementary optical and near-IR high angular resolution HST images have also been used. Compared to classical optical long-slit 
analysis, our data provides information about the entire system, and not only in pre-selected orientations. 

Results. The extinction exhibits a very complex and patchy structure in ULIRGs on kpc scales, from basically transparent regions to 
others deeply embedded in dust (Av)-O.O to Av-8.0 mag). Nuclear extinction covers a broad range in Ay from 0.6 to 6 mag, 69% of 
the nuclei having Av>2.0 mag. Extinction in the external regions is substantially lower than in the nuclei with 64% of the ULIRGs 
in the sample having median Ay of less than 2 mag for the entire galaxy. While post-coalescence nuclei tend to cluster around Ay 
values of 2 to 3 mag, pre-coalescence nuclei appear more homogeneously distributed over the entire 0.4 mag <Ay < 7.7 mag range. 
For the average extinction (Ay^2.0) derived for the ULIRGs of the sample, the ratio of the de-reddened to observed SFR values is 6. 
The extinction-corrected, Hor-based SFR ranges from 10 to 300 MO yr~^ For only 28% of the cases the de-reddened SFR is <20 MO 
yr~\ whereas for the observed SFR this percentage increases to 72%. The IR-based SFR is always higher than the optical-based one, 
with differences ranging from about 2 to up to 30. The nuclear observed SFR has an average contribution to the total one of 16% for 
the entire sample. Once corrected for extinction, the average value becomes 31%. 

Because of mostly extinction effects, the optical (I-band) half-light radius in the sample galaxies is on average a factor 2.3 larger than 
the corresponding near-IR (H-band) value. 
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1. Introduction 

The significance of galaxies with an infrared excess was high- 
lighted early by Low & Kleinmann (1968) and Rieke & Low 
(1972), and lately recognized by observations of the IRAS satel- 
lite (Soifer et al. 1984). Among the sources that emit the bulk 
of their energy in the infrared, the Ultra-Luminous Infrared 
Galaxies (ULIRGs, 10^^Lq< Lboi~ Lir[8 - lOOOyum] < lO^^L©), 
are amongst the brightest objects in the local Universe (see re- 
views in Lonsdale et al. 2006; Sanders & Mirabel 1996). The 
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process of merging with accompanying starbursts plays a key 
role in these galaxies, and appears to transform spirals into low- 
intermediate mass ellipticals (e.g., Colina et al. 2001; Genzel et 
al.2001; Tacconi et al.2002; Naab et al.2006 and references 
therein). This star formation activity is the main energy source in 
ULIRGs, although in some cases the contribution of an obscured 
AGN is relevant (see Nardini et al.2008; Risaliti et al.2006). 
These galaxies contain large quantities of gas and dust that, inde- 
pendent of the nature of the energy source, are heated re-emitting 
most of their energy in the infrared. 

High angular resolution HST optical and near-infrared im- 
ages (Farrah et al. 2001 ; Bushouse et al. 2002) infer that ULIRGs 
have in general a complex stellar structure. On the one hand, 
these structures are caused by the merger itself that transforms 
the stellar distribution in the parent galaxies (i.e., into tidal tails, 
nuclear and extended star-forming regions, and double nuclei) 
as they evolve through the different phases of the interaction. 
On the other hand, large amounts of gas and dust settle into ex- 
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tended lanes and/or filaments creating additional structures due 
to non-uniform extinction eff'ects (see V, I, and H band images in 
Farrah et al. 2001 and Bushouse et al. 2002). The impact of the 
dust eff'ects is lower in the near-IR, allowing us to identify the 
true stellar morphology of the galaxies (e.g., HST observations 
from Scoville et al. 2000 and Bushouse et al. 2002). This is em- 
phasized by studies based on HST color maps (F814W-F606W), 
which detect a non uniform color structure with blue and red 
compact knots that are most likely caused by unobscured star 
formation, and dust-enshrouded star formation or remnant nu- 
clei from the progenitor galaxies (Surace et al. 1998; Farrah 
et al. 2001). Similar results have been obtained for luminous 
infrared galaxies (LIRGs, 1O^^L0< Lboi~ Lir[8 - lOOOyum] < 
10 ^^Lq) using near-IR HST images. Complex dust features in 
their centers, with an average optical extinction (Ay) of between 
3 and 5 mag, are also observed (Alonso-Herrero et al. 2006). 

All the afore mentioned HST imaging studies indicate that 
although the dust tends to be concentrated in the inner few kpc, 
the global distribution is very patchy on scales of hundreds to 
thousands of parsecs, and therefore non-uniform extinction ef- 
fects play a major role. ULIRGs have been intensively studied 
with narrow, long-slit spectroscopy (e.g., Veilleux et al. 1995), 
but despite its importance, the spectroscopic derivation of the 
two-dimensional dust and gas distribution has not been the sub- 
ject of much scrutiny. Previous integral field spectroscopy (IFS) 
studies of specific targets demonstrated the lack of coincidence 
between the true, dynamical nucleus of the galaxy and the re- 
gion identified as the optical nucleus, which can be separated 
by distances from 0.5 (Arp 299, Garcia-Marm et al. 2006) to 1.3 
kpc (IRAS 17208-0014, Arribas & Colina 2003) . Other studies 
have also shown the large diff'erences (of up to a factor of 10) in 
extinction between the nuclear, circum-nuclear, and extranuclear 
regions due to the patchy and non-uniform distribution of the 
dust on kpc-scales (e.g. IRAS 121 12-0305, Colina et al. 2000; 
Arp 299, Garcia-Marm et al. 2006). These results emphasize the 
importance of knowing (and correcting for) the two-dimensional 
structure of the extinction before using rest-frame optical lines 
(mainly Ha) or structural parameters (e.g., eff'ective radius) in 
the derivation of star formation rates and dynamical masses in 
these galaxies, and their intermediate- and high-z analogs. 

This is the second paper, in a series aimed at studying in 
detail the internal structure and kinematics of local ULIRGs 
(Garcia-Marm et al. 2009, Paper I). Here, we present the first 
systematic two-dimensional analysis of the dust extinction in a 
representative sample of low-z ULIRGs using optical integral 
field spectroscopy complemented with existing archival HST 
high angular resolution images. The paper is organized as fol- 
lows. The galaxy sample is presented in Sect. 2, whereas Sect. 

3 is dedicated to the observations and data reduction, and Sect. 

4 presents the data analysis. Section 5 discusses the extinction 
structure of ULIRGs on kpc scales, and the implications for the 
derivation of star formation rates and dynamical masses in low- 
and high-z ULIRGs. Finally, in Sect. 6, a brief summary of the 
main results is given. 

The current analysis that we are carrying out in this sam- 
ple of ULIRGs is part of a larger survey that is investigating 
the two-dimensional extinction, ionization, and kinematic kpc- 
scale structure of a representative sample of low-z LIRGs and 
ULIRGs, and its implications for their high-redshift analogs. The 
survey is based on the use of optical IFS data obtained with 
diff'erent facilities. The study of northern hemisphere ULIRGs 
has been performed mostly with the integral field unit (IFU) 
INTEGRAL (Arribas et al. 1998). Several papers presenting the 
results for individual (Colina, Arribas & Borne, 1999; Garcia- 



Marm et al. 2006) or small subsamples of galaxies (see Colina et 
al. 2005; Monreal-Ibero et al. 2007 and references therein) have 
already been published. 

To include southern systems and extend the sample to 
lower luminosity objects covering most of the (U)LIRGs range 
(log(LiR/Lo) = 11.0-12.6), we observed an additional sample of 
about 40 systems using VIMOS IFU (LeFevre et al. 2003). The 
VIMOS sample description and first results were presented in 
Arribas et al. (2008). 

We also observed the majority of the northern hemi- 
sphere galaxies from the volume-limited (distances of about 35- 
75 Mpc) local sample of LIRGs presented in Alonso-Herrero et 
al. (2006); further details of the selection criteria can be found 
in Alonso-Herrero et al. 2006. These galaxies, with an average 
luminosity of log(LiR/Lo)= 11.32, were observed with the IFU 
PMAS (Roth et al. 2005). The data and first results are presented 
in Alonso-Herrero et al. (2009, Paper I.). 

In addition to this, near-IR IFS with SINFONI of a subsam- 
ple of LIRGs and ULIRGs is also being conducted (see Bedregal 
et al. 2009 for first results) to investigate the multi wavelength 
(optical and near-IR) properties of these galaxies. 

Throughout this paper, we use Da=0.7, Dm=0.3, and //o=70 
km s"^ Mpc"^ 

2. The sample of ULIRGs 

Our original sample of ULIRGs with available optical 
INTEGRAL IFS consists of 22 systems selected to be repre- 
sentative of this galaxy class (details on the sample selection 
can be found in Garcia-Marm et al. 2009, Paper I). These sys- 
tems exhibit a large variety of morphologies, from wide pairs to 
close pairs with well developed tidal tails to single nucleus merg- 
ers. The sample covers the luminosity range 11.8<log(LiR/Lo) 
<12.6, and includes all the diff'erent classes of nuclear activity 
(Hii-, LINER-, and Seyfert-like). A few galaxies of the sample 
are not present in the IFS-based analysis of the dust distribution 
for a variety of reasons, such as the lack of detection or non cov- 
erage of HyS (Arp 220 and IRAS 09427 -hl929), low S/N, which 
precludes the detection of emission lines (IRAS 13469+5833), 
limited spatial resolution (IRAS 13342+3932), and severe AGN 
contamination (Mrk 231). Using the projected distance between 
the nuclei of the parent galaxies, the sample has been further 
divided into two categories, pre- (nuclear distance >1.5 kpc) 
and post-coalescence (nuclear distance <1.5 kpc). This separa- 
tion allows us to investigate and discriminate the characteristics 
of galaxies in the final phases of the merger from those in ear- 
lier dynamical stages. Given this, the two-dimensional dust dis- 
tribution has been derived for 17 (U)LIRGs, nine of them pre- 
coalescence systems and eight post-coalescence (Table [B. 

3. Observations and data reduction 

IFS data of the galaxy sample was obtained between 1998 and 
2004 using INTEGRAL, a fiber-based optical integral field sys- 
tem (Arribas et al. 1998) connected to the Wide Field Fibre Optic 
Spectrograph (WYFFOS; Bingham et al. 1994) and mounted on 
the 4.2 m William Herschel Telescope. Depending on the struc- 
ture and compactness of the ULIRGs,three diff'erent INTEGRAL 
configurations were used: the so-called standard bundles 1, 2, 
and 3 (SBl with fiber diameter 0.45 arcsec and field-of-view 
(FoV) 7.8x 6.4 arcsec^, SB2 with fiber diameter 0.9 arcsec and 
FoV 16. Ox 12.3 arcsec^, and SB3 with fiber diameter 2.7 arc- 
sec and FoV 33. 6x 29.4 arcsec^). The spectra were taken with 
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^ The galaxies are listed according to decreasing nuclear separation (see Garcia-Marm et al 2009). 
^ Lir(8-1000 \im) was derived following Sanders & Mirabel 1996) 

^ Median extinction of the entire system with the deviation that effectively represents the distribution of the extinction values over the entire area 
where the Hof and Hy5 measurements can be obtained. For pre-coalescence systems with one pointing per nuclei, the median has been calculated 
for each pointing. 

^ Nuclear extinction of the galaxy. The uncertainties are about 20% 

^ Maximum extinction of the system. 

^ Separation between the optical nucleus, identified as the peak position in the F814W band of the HST images, and the position of the peak of 
extinction. 

^ Eff'ective radius measured from the HST WFPC2 camera, with the filter F814W, corresponding to a central wavelength of 0.82 yum and a filter 
width of 0.17yL/m. 

'^ Eff'ective radius measured from the HST NICM0S2 camera, with the filter F160W, corresponding to a central wavelength of 1.55 \im and a 
filter width of 0.40 yt/m. 
' IP means interacting pair, DN double nucleus and SN single nucleus. 

^ Eff'ective radius measured from the HST WFPC2 camera, with the filter F606W, corresponding to a central wavelength of 0.58 yum, and a filter 
width of 0.16yL/m. 

^ Eff'ective radius measured from the HST WFPC2 camera, with the filter F439W, corresponding to a central wavelength of 0.43 yt/m and a filter 
width of 0.05 yL/m. 

' Eff'ective radius measured from the HST ACS camera, with the filter F435W, corresponding to a central wavelength of 0.43 yum and a filter 
width of 0.14 yL/m. 



a 600 lines mm ^ grating, providing an eflTective spectral res- 
olution (EWHM) of approximately 6.0, 6.0 and 9.8 A for the 
SBl, SB2, and SB3 bundles, respectivel}Q. The covered spectral 
range of interest was /1/14500-7000 A rest- frame. The reduction 
and calibration of the lES data were performed inside the IRAlQ 

^ These values correspond to the old camera mounted on WYFFOS, 
that was used for the present observations. From August 2004 a new 
camera was commissioned for the instrument. See more details in 
http : //w w w. iac . es/proy ecto/integral 

^ The IRAF software is distributed by the National Optical 
Astronomy Observatory (NOAO), which is operated by the Association 
of Universities for Research in Astronomy (AURA), Inc., in coopera- 
tion with the National Science Foundation. 



environment, and followed the standard procedures applied to 
this type of data (see Arribas et al. 1997 and references therein). 
Full details of the observations and data reduction can be found 
in Garcia-Marin et al. (2009). 

Complementary HST archive images were also used. 
Specifically, we used optical WEPC2 E439W, E606W, and 
E814W data, available for 12, 12, and 100% of the ULIRGs un- 
der study, and near-IR NICMOS E160W data, available for 60% 
of the sample. Filters E439W, F814W, and E160W, are the HST 
analogs to the ground-based Johnson-Cousins B, I, and H, re- 
spectively (Origlia & Leitherer 2000). The filter F606W is equiv- 
alent to a wide V- filter. In addition, E435W ACS images, equiv- 
alent to the Johnson-B filter and available for 29% of the sam- 
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pie, were also used. The HST images were calibrated on the fly, 
with the best available reference files. In the case of NICMOS 
data, the combination of the dithered individual exposures was 
repeated. All images were sky (or background) subtracted. 

4. Data analysis 

In all galaxies with available information, the two-dimensional 
structure of the extinction/dust was derived using the Hof/HyS line 
ratio (B aimer decrement), a foreground dust screen model and 
a mean interstellar extinction law based on Savage & Mathis 
(1979). The line fluxes were obtained by fitting each emission 
line to a single Gaussian function using the DIPSO package 
(Howarth & Murray 1988) inside the STARLINK environment 
The presence of underlying stellar hydrogen absorption lines 
was not detected and therefore no correction was applied when 
measuring the flux of the corresponding emission lines. This 
would introduce a small overestimate in the values derived for 
the optical extinction since the equivalent widths of the hydro- 
gen emission lines are larger than 40 A in the extended high 
surface brightness regions where they are detected, much larger 
than the EW^\)^(Y{a)=EW^\)^(Y{p)--2 A of the underlying stars. 

It is known that difl'erent emission line ratios infer difl'er- 
ent extinction values, indicating that the observed recombina- 
tion lines originate at difl'erent depths in these dusty regions. For 
this reason, especially in dust-enshrouded environments such as 
the central regions of ULIRGs where large quantities of dust 
are detected, the B aimer decrement provides only a lower limit 
to the true extinction. A clear example of this is the nucleus 
of the galaxy IC 694, a member of the pre-coalescence system 
Arp 299. The nuclear extinction derived using the B aimer decre- 
ment is Ay -3.0, whereas with the Paa/Ha line ratio it is Ay -6.0 
(Garcia-Marm et al. 2006). 

The HST ACS, WFPC2, and NICMOS images were used to 
obtain the optical (B or V and I bands) and near-IR (H band) 
eff'ective radius for each of the sample galaxies, defined as the 
radius that encloses half the luminosity of the galaxy. It was de- 
rived using sky- or background-corrected images, using aper- 
tures of increasingly large radius being centered on the optical 
nucleus of each galaxy. For the pre-coalescence systems, con- 
tamination from neighboring galaxy was prevented by the use of 
masks. In all cases, we include the outer low surface brightness 
regions of the galaxies, until the signal was at the same level as 
the background. 

5. Results and discussion 

5.1. Two-dimensional Ha/Hj3-based extinction in pre- and 
post-coalescence systems 

The two-dimensional Ha/HjS-bsLScd extinction maps of the 
sample galaxies cover areas of several kpc on a side, rang- 
ing from about 2.5 kpc in highly obscured systems (e.g., 
IRAS 17208-0014) to 10-12 kpc in galaxy pairs (e.g., IRAS 
08572-h3915 or IRAS 14348-1447). The angular resolution of 
the data provides extinction values on scales of between 0.7 and 
3 kpc, depending on the redshift of the source (see specific val- 
ues for each galaxy in Table[T]). On these scales, the extinction 
maps (Fig. [T]) show that the dust in ULIRGs is not uniformly 
distributed. It shows instead a patchy structure including almost 
transparent regions with very low dust content, and other regions 

^ S ee |http : //www. starlink.rl . ac .u 



that are highly absorbed where, in the optical domain, only a 
lower limit to the extinction can be derived. 

The visual extinction (Av=3.1xE(B-V)) ranges from 
Av<0.2 mag, generally measured in the external regions, to 
Ay -9.0 mag, measured in the central kpc of of IRAS 17208- 
0014. The highest extinctions are usually being measured in the 
nuclear regions, whose values cover a broad range of Ay be- 
tween 0.6 and 7.7 mag, 69% of the nuclei having Av>2.0 mag. 
These high extinction values are similar to those measured in the 
Paof emitting regions of LIRGs, whose average extinctions are 
found to be Ay =3 - 6 mag (Alonso-Herrero et al. 2006) but lower 
than the values derived for ULIRGs using the near-IR Pao^/Bry 
line ratio, where visual extinctions (Ay) in excess of 10 magni- 
tudes are relatively common in ULIRGs (Murphy et al. 2001). 
Both optical and near-IR extinction values are also systemati- 
cally lower than the Spitzer mid-IR based ones derived for local 
ULIRGs, which are found to be <20 - 30 mag for PAH-emitting 
normal starbursts, and »30 mag for centrally concentrated en- 
ergy sources with no PAH emission (Imanishi et al. 2007). 

Pre- and post-coalescence systems are in diff'erent phases of 
the evolutionary merging process, and therefore their dust and 
gas content and distribution could reflect this. Pre- and post- 
coalescence systems do have similar nuclear extinctions with 
average values of 3.0 and 2.5 mag, respectively. However, while 
while extinction values for post-coalescence nuclei tend to be 
between 2 and 3 mag of visual extinction, the values for pre- 
coalescence nuclei appear to be more widely distributed over the 
entire 0.5 mag < Ay < 8 mag range (Fig. [2] center and Table 1 for 
specific values). 

The median extinction derived for entire galaxies (Fig. [2] 
right) show a clear indication that extinction in the circumnu- 
clear and external regions is substantially lower than in the nu- 
clei. While only 31% of all nuclei have visual extinctions of 2 
mag or less, 60% of the galaxies do have median extinctions with 
those low values. Of all the systems, pre-coalescence ULIRGs 
show on average a diff'erence of 1.2 mag between the median 
visual extinction for the entire galaxy (1.8 mag) and the nu- 
clear regions (3.0 mag). Post-coalescence systems do however 
show similar values (2.1 and 2.5 mag for median and nuclear 
extinctions, respectively). This behavior is clearly explained as 
a consequence of the diff'erent morphology of pre- and post- 
coalescence systems. Pre-coalescence systems contain several 
low-extinction Ha high-surface brightness regions at distances 
of several kpc away from the nucleus associated with star for- 
mation along the tidal tails and candidates to tidal dwarf galaxies 
(Monreal-Ibero et al. 2006). Post-coalescence systems are on the 
other hand more compact and dominated by the emission from 
the nucleus. The importance of these eff'ects in deriving star for- 
mation rates in dust-enshrouded systems is discussed next. 

5.2. Implications of the extinction structure in the derivation 
of the Star Formation Rates 

Star formation rates (SFR) in galaxies are derived using diff'erent 
tracers, from ultraviolet light and optical forbidden lines ([On] 
3727A), to hydrogen recombination lines (e.g. Ha, PacK, Bry), 
and IR broad-band flux (see e.g. Kennicutt 1998). The Ha line 
is the most intense emission line in the optical, and shifts into 
the near- and mid-IR spectral range for intermediate and high-z 
galaxies. Because of that, it is one of the most common spec- 
tral features used in deriving SFR in all kinds of galaxies at all 
redshifts. 
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For obscured and dust-enshrouded galaxies such the low-z 
ULIRGs and the high-z submillimeter and 24 jim bright Spitzer 
galaxies, the extinction correction applied has important conse- 
quences in the reliability of the Ha-based SFR. The SFR and Ha 
luminosity are directly proportional as given by the expression 

SFR(M^yr-') = 7.9 x IQ-'^Luaierg s'') (1) 

for a continuous star formation, a Salpeter initial mass function 
with upper and lower mass limits of 100 and 0.1 M©, respec- 
tively, and assuming that all potentially ionizing photons ionize 
the surrounding interstellar medium (Kennicutt 1998). Most dust 
extinction in low-z ULIRGs to date are obtained through narrow 
(typically 1 to 1.5 arcsec) long-slit spectra along a given orienta- 
tion centered on the nucleus of the galaxy (if post-coalescence), 
or connecting the nuclei of interacting pairs (if pre-coalescence). 
Observational difficulties arising from differential atmospheric 
refraction and positioning of the galaxy on the slit translate into 
an increased uncertainty in the measured values. Since the width 
of the slit does not cover the entire Ha emitting regions (see 
Figs. 1-4 for angular extent and structures), the extinction val- 
ues derived from these data will be dominated by the emission 
from the nucleus. This would underestimate the Ha luminosity, 
and therefore the SFR, by an average factor of about 4 and 7 in 
the pre- and post-coalescence ULIRGs, respectively (see Table 
2 for the specific ratios of the observed nuclear to total Ha flux 
(Col. 3), and the corresponding extinction corrected luminosities 
(Cols. 5 and 8) for each galaxy). 

For the typical median ULIRG extinction (Ay^2.0 see Fig. [2] 
right panel), the use of observed Ha luminosities underestimates 
the SFR by a factor 6 with respect to the extinction-corrected 
value. Only 28% of the galaxies in the sample have de-reddened 
SFR of less than 20 Mo yr"\ whereas for the uncorrected SFR 
this percentage increases to 72%. The total de-reddened SFR 
of the present ULIRGs sample ranges from about 10 to 300 
Mo yr"\ for a similar range of SFR values for pre- and post- 
coalescence systems. Finally, we note that because of the lower 
S /N ratio of the HJ3 line, the extinction maps do not cover the 
same field as the Ha maps (see from Figs. 1 to 4). This would 
lead to an additional underestimate of the two-dimensional ex- 
tinction corrected Ha fluxes (and thus of the SFR values). In any 
case, and since this mainly afl'ects the outer low-surface bright- 
ness regions, no significant changes should be expected. 

The SFR is also found to be proportional to the Ljr. 
Considering continuous bursts of 10-100 Myr and applying the 
same range of IMF parameters as used in Eq. (1) above, the pro- 
portionality is given by Kennicutt (1998) 

SFRiM^yr'^) = 4.5 x lO-"^"^ Li^ierg s'^). (2) 

The SFR as derived from the IR luminosity is systematically 
higher than that obtained from the Ha measurements (see Table 
3). The ratio S FR(IR) / S FR(L(H a) con) covers a wide range, 
from about 2 to 10, and tends to be smaller in galaxies with 
higher nuclear extinctions, i.e., IRAS 121 12-h0305, IRAS14348- 
1447, and IRAS 1 6007 -h3743. This is because of the large 
amounts of dust obscuring the nuclear regions, which cause the 
innermost regions to be more optically thick to optical than far 
infrared radiation. Thus, even by applying the two-dimensional 
extinction corrections to the observed Ha fluxes, the derived 
SFRs would represent 10% to 50% of the IR-based SFR value, 
which is assumed to be free of extinction efl'ects, and therefore 
represent the true SFRs. 



5.3. Implications of the extinction structure for the derivation 
of the dynamical mass 

Dynamical masses (oc cr^ x 7?hm) of galaxies are often derived 
from measurements of the velocity dispersion of the stars or gas, 
and assuming the radius (Rhm) where half of the mass is con- 
tained, is known. While the velocity dispersion can be directly 
obtained, to first order, from the full-width-half-maximum of the 
absorption/emission line profiles, the half-mass radius is not a di- 
rect observable and has to be inferred from the half-light radius. 

The most common lines used to derive velocity dispersions 
in low-z ULIRGs and their high-z analogs are in the optical (e.g.. 
Hp and Ha) and near-IR (e.g., Bry, CO stellar absorption bands) 
rest- frame, and are therefore aff'ected by non-uniform extinction 
(see Sect. 5.1). The behavior of the extinction in ULIRGs with 
large values in the nuclear regions and an outward decreasing 
gradient in the circumnuclear and extranuclear regions (see Sect. 
5.1 and Fig. 1), produces an apparent stellar light distribution 
that diff'ers from that of the true mass distribution. On the one 
hand, for a given rest- frame wavelength, the half-light radius will 
always be an upper limit to the true R^m as the light from stars in 
the circum- and extranuclear regions are less obscured than those 
in the nuclear regions. On the other hand, the extinction has a 
strong wavelength dependence becoming higher towards bluer 
wavelengths. For a given galaxy, the half-light radius derived in 
the optical will always be systematically larger than those de- 
rived in the near-IR. This eff'ect is observed in our sample galax- 
ies where for those galaxies with available HS T optical (0.8//m, 
WFPC2 F814W) and near-IR (1.6 yum, NICM0S2 F160W) im- 
ages (see Table 1), the optical half-light radius is on average 2.2 
times the corresponding value in the near-IR, in some systems 
up to 5 times larger. Therefore, in obscured systems such as lo- 
cal and high-z ULIRGs, the use of optical (red) rest- frame light 
will produce an average overestimate of the dynamical mass by 
a factor of 2.2 merely because of extinction eff'ects. This factor 
could be even larger if eff'ective radii are measured in the blue 
(B) or visual (V) rest- frames (see specific values in Table 1), as 
frequently performed for in high-z galaxies. 

In addition to these extinction eff'ects, other eff'ects such as 
stellar population gradients can also be important. Since young 
(10 Myr or less) massive stars and older (100 Myr to Gyrs) stel- 
lar populations have mass-to-light ratios that diff'er by up to two 
orders of magnitude in the optical and near-IR (e.g. Starburst99, 
Leitherer et al. 1999), their non-uniform spatial distribution is 
also relevant to evaluating the uncertainty in the measured half- 
light radius, and therefore in the dynamical mass. With only two 
filters available for most of the sample galaxies, a detailed anal- 
ysis of these eff'ects is beyond the scope of the present work. 



6. Summary 

This paper has presented the study of the two-dimensional ex- 
tinction structure of a representative sample of 17 ULIRGs, 9 
pre-coalescence (i.e., interacting pairs with projected nuclear 
separations of at least 1.5 kpc), and 8 post-coalescence (either 
single nucleus or double nuclei with projected separations of 
less than 1.5 kpc). The extinction maps are based on the mea- 
surement of the Ha/H/3 line ratio obtained from integral field 
spectroscopy with the instrument INTEGRAL on the William 
Herschel telescope. In comparison with more classical narrow 
long-slit spectroscopy, this study infers the two-dimensional ex- 
tinction structure of ULIRGs on kpc scales over areas of a few to 
several kpc on a side. The main results are summarized below: 
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Table 2. Observed Ha fluxes and extinction corrected luminosities for the sample galaxies. 



Galaxy 


^(HQf)ob^ 


-/^(Hof ^j^uclear /total 


^(HQf)ob^ 


-^(^-H cxtcor^ 








IRAS 13156+0435N 


0.9 


0.12 


3.0 


21.0 


21.1 


0.36 


4.6 


IRAS 13156+0435S 


1.2 


0.09 


3.9 


9.0 


6.1 


0.35 


1.3 


IRAS 18580+6527 


2.5 


[E] 0.01 


21.6 


78.0 


68.2 


0.22 


0.4 






[W] 0.28 








6.05 


15.9 


IRAS 16007+3743 


2.6 


[E] 0.03 


25.1 


152.0 


134.7 


0.75 


17.7 






[W] 0.22 








5.52 


97.0 


IRAS 06268+3509 


0.6 


[N] 0.09 


4.7 


121.0 


28.0 


0.42 


66.5 






[S] 0.18 








0.85 


4.4 


IRAS 08572+3915 


5.8 


[N] 0.08 


4.6 


15.5 


21.1 


0.37 


2.0 






[S] 0.05 








0.23 


0.3 


IRAS 14348-1447 


15.0 


[N] 0.03 


25.6 


381.0 


333.0 


0.77 


52.3 






[S] 0.08 








2.05 


31.5 


Mrk 463 


36.9 


[E] 0.25 


21.8 


56.0 


48.2 


5.45 


10.4 






[W] 0.02 








0.44 


1.8 


Arp 299/NGC 3690 


181.0 


0.03^;, 


4.0 


29.0 


21.5 


0.12 


0.9 


Arp 299/IC 694 


89.0 


0.03" 


2.0 


16.0 


11.9 


0.06 


6.2 


IRAS 12112+0305 


5.8 


[N] 0.08 


7.6 


213.0 


34.0 


0.60 


9.0 






[S] 0.02 








0.15 


155.0 


IRAS 06487+2208 


2.5 


0.19 


13.9 


117.0 


96.3 


2.64 


23.3 


IRAS 11087+5351 


0.8 


0.12 


4.4 


11.0 


12.0 


0.53 


7.4 


Mrk 273^ 




0.08 












IRAS 12490-1009 


1.0 


0.13 


2.6 


12.0 


8.1 


0.34 


0.9 


IRAS 14060+2919 


2.3 


0.18 


8.2 


45.0 


30.8 


1.48 


10.7 


IRAS 15206+3342 


42.7 


0.38 


171.8 


322.0 


348.7 


65.3 


158.5 


IRAS 15250+3609 


4.4 


0.26 


3.2 


87.5 


25.0 


0.83 


4.4 


IRAS 17208-0014 


17.9 


0.10 


7.7 


160.2 


j 


0.7 


43.0 



^ Integrated Ha observed flux, with units x 10 erg s ^ cm ^. Note that in the cases of IRAS 13 156+0435 and Arp 299 we derived the integrated 
flux for each galaxy of the interacting pair using their individual pointings. 
^ Ratio between the Ha nuclear, identified as the peak position in the optical continuum and total observed fluxes. 
^ Observed Ha luminosity derived from the integrated flux in column 2. Units in xlO"^^ erg s~^ 
^ Integrated two-dimensional extinction corrected Ha luminosity, with units in xlO^^ erg s~^ 

^ Integrated Ha extinction-corrected luminosity, derived using the total observed flux and the median extinction, with units in xlO"^^ erg s~^ 
^ Observed nuclear Ha luminosity with units in xlO"^^ erg s~^ 

^ Extinction corrected nuclear Ha luminosity, derived using the nuclear extinction of each individual galaxy. Units in xlO"^^ erg s~^ 

^ For this particular galaxy we have calculated the ratio using the position of the true nucleus of the galaxy, which does not coincide with the 

optical one (Garcia-Marm et al 2006). 
' The data for Mrk 273 do not have an absolute flux calibration. Hence, only the relative measurement between the nuclear and integrated Ha 

flux is given here. 

^ As shown in Fig. lb, the area of this galaxy with available extinction values is limited, also presenting a vey large peak value. The use of this 
value for correcting the integrated flux is not valid and therefore is not included in the table. 



- In agreement with previous studies, the analysis of our data 
has detected a very complex and patchy extinction struc- 
ture in ULIRGs on scales of kpc, from basically transpar- 
ent regions to others so deeply embedded in dust that only a 
lower limit to their extinction can be derived (Ay -0.0 mag 
to Ay -8.0 mag). Nuclear extinction values of the present 
galaxy sample cover a broad range in Ay from 0.6 to 8 mag, 
69% of the nuclei having Ay>2.0 mag. Extinction in the ex- 
ternal regions is substantially lower than in the nuclei with 
60% of the ULIRGs in the sample having a median Ay for 
the entire galaxy of less than two magnitudes. 

- ULIRGs classified as pre- and post-coalescence exhibit sig- 
nificant diflTerences in their extinction properties. While post- 
coalescence nuclei tend to cluster around Ay values of 2 to 
3 mag, pre-coalescence nuclei appear more homogeneously 
distributed over the entire 0.5 mag < Ay < 6 mag range. 
Pre-coalescence ULIRGs show on average a diflTerence of 
1.2 mag between the nuclear (Ay of 3.0 mag) and me- 
dian extinction for the entire system (Ay of 1.8 mag), while 
post-coalescence systems do show similar values (2.1 and 



2.5 mag for median and nuclear extinctions, respectively). 
This behavior is explained as a consequence of the diflfer- 
ent distribution of the star-forming regions in pre- and post- 
coalescence systems. Pre-coalescence systems contain sev- 
eral low-extinction Ha high-surface brightness regions at 
distances of several kpc away from the nucleus associated 
with star formation along the tidal tails. Post-coalescence 
systems are more compact and dominated by the emission 
from the nuclear regions. 

- Considering the median extinctions measured in the sample 
galaxies (Ay^2.0), the SFR values based on the de-reddened 
Ha fluxes increase by a factor 6 with respect to those de- 
rived from the observed Ha fluxes. The HcK-based extinction- 
corrected SFR ranges from about 10 to 300 Mo yr"^ . The de- 
reddened SFR is less than 20 Mo yr"^ in 28% of the galaxies, 
whereas for the observed SFR this percentage increases up to 
72%. No significant diflferences in the SFR have been found 
between the pre- and post-coalescence systems. 

- In all galaxies, the IR-based SFR is higher than the de- 
reddened Ha value. The ratio ranges from 2 to 10, with 
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Galaxy logiLi^L SFR(L{UaU,t S FR(L{Ua)^,,_,,Q S FR(L{Ua\,,,t S FRjlR^ 



IRAS 13156+0435 


12.13 


5.5 


26.2 


23.7 


232.9 


IRAS 18580+6527£ 


12.26 


17.1 


53.8 


61.6 


318.7 


IRAS 16007+3743" 


12.11 


19.8 


106.4 


120. 


222.5 


IRAS 06268+3509 


12.51 


3.7 


22.1 


95.6 


493.5 


IRAS 08572+3915 


12.17 


3.6 


16.7 


12.2 


255.5 


IRAS 14348-1447 


12.39 


20.2 


263.1 


301.0 


424.1 


Mrk 463 


11.81 


17.2 


38.0 


44.2 


111.5 


Arp 299 


11.81 


4.7 


26.3 


35.3 


111.5 


IRAS 12112+0305 


12.37 


6.0 


26.9 


168.3 


405.0 


IRAS 06487+2208 


12.57 


11.0 


76.0 


92.4 


641.8 


IRAS 11087+5351 


12.13 


3.5 


9.5 


8.7 


233.0 


Mrk 273^ 


12.18 








216.5 


IRAS 12490-1009 


12.07 


2.0 


6.4 


9.5 


203.0 


IRAS 14060+2919 


12.18 


6.4 


24.3 


35.5 


261.5 


IRAS 15206+3342 


12.27 


135.7 


275.5 


254.3 


321.6 


IRAS 15250+3609 


12.09 


2.5 


19.7 


69.1 


212.5 


IRAS 17208-0014 


12.43 


6.1 




126.5 


465.0 



^ IR luminosity derived following Sanders & Mirabel (1996) 

^ SFR derived using the total observed Ha flux given in column 2 of Table 2. Units are MO yr~^ 

^ SFR derived using the Ha median extinction corrected values as specified in Table 2 column 6. Units are Mo yr~^ 

^ SFR calculated using the two-dimensional extinction corrected Ha luminosity as given in Table 2 column 5. Units are MO yr~^ 

^ SFR based on the IR-based luminosity values as given in column 2 of the present table. Units are Mo yr~^ 

^ Strong evidence for a Seyfert 2 nucleus, for which we have not corrected this values 



the smaller differences in galaxies with higher extinction in 
their nuclear regions. Thus, assuming IR luminosities mea- 
sure the true SFR in these systems, SFRs based on full two- 
dimensional extinction corrected Ha luminosities recover 
between 10% and 50% of the true SFR. 
- Because of the high nuclear extinctions and the outward de- 
creasing extinction gradient, the optical (I-band) half-light 
radius in the sample of ULIRGs is on average a factor 
2.2 larger than the corresponding value in the near-IR (H- 
band), and in some systems even a factor 5 larger. Dynamical 
masses in ULIRGs and their high-z analogs would be on 
average overestimated by a factor 2.2 due merely to extinc- 
tion effects if the effective radius of the galaxy is derived di- 
rectly from observed rest- frame optical images. This factor 
could be even larger if bluer broad-band filters (i.e., V and 
B bands) were used to measure the effective radius as often 
the case with high-z galaxies. Additional stellar population 
effects could also be relevant because of their different mass- 
to-light ratios, as should be investigated by future detailed 
studies. 
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Fig.l. a. From left to right, HST WFPC2/F814W image, Ua and extinction Ay INTEGRAL maps for a set of ULIRGs. The 
extinction maps were derived using the Ha/HyS line ratio. The HST and Ha maps are represented in logarithmic scale, whereas for 
the extinction map the scale is linear. As a reference, the white contours shown on every map represent different Ha intensity levels 
of each galaxy. For the Ha map, we included a color code, which is given in relative flux units. Each pixel of the INTEGRAL maps 
represents the fiber size of the fiber used ((y.'9). The horizontal scale represents 5 kpc. Galaxies orientation is north up, east to the 
left. 
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Fig. 1. b. Same as Fig.[T^, but for a different set of ULIRGs. 
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Fig. 1. c. Same as Fig.[T^, but for Arp 299, an interacting galaxy formed by NGC 3690 (West) and IC 694 (East). In this case the 
scale represents 1 kpc. 




Aa(arcsec) Aa(arcsec) Aa(arcsec) 

Fig. 1. d. Same as Fig.[T^, but for two different ULIRGs observed with a fiber size of Off 45. 
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Fig. 2. Distribution of the extinction values in the sample of ULIRGs. The black dashed histogram indicates pre-coalescence sys- 
tems. The red dashed histogram indicates post-coalescence galaxies. Left panel: Distribution of the sample of galaxies according 
to the peak of extinction for each ULIRG. Center panel: Galaxy distribution according to the extinction in the nuclear regions. 
Right panel: Distribution of the sample of galaxies according to the median extinction of the systems. Note that in the cases of 
IRAS 13156-h0435 and Arp 299, we derived the median for each galaxy of the interacting system using the individual pointings 
with the SB2 bundle. 



